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Abstract — Traditional radar sensing typically involves matched 
filtering between the received signal and the shape of the 
transmitted pulse. Under the confinement of classic sampling 
theorem this requires that the received signals must first be 
sampled at twice the baseband bandwidth, in order to avoid 
aliasing. The growing demands for target distinction capability 
and spatial resolution imply significant growth in the bandwidth 
of the transmitted pulse. Thus, correlation based radar systems 
require high sampling rates, and with the large amounts of 
data sampled also necessitate vast memory capacity. In addition, 
real-time processing of the data typically results in high power 
consumption. Recently, new approaches for radar sensing and 
detection were introduced, based on the Finite Rate of Innovation 
and Xampling frameworks. These techniques allow significant 
reduction in sampling rate, implying potential power savings, 
while maintaining the system's detection capabilities at high 
enough SNR. Here we present for the first time a design and 
implementation of a Xampling-based hardware prototype that 
allows sampling of radar signals at rates much lower than 
Nyquist. We demostrate by real-time analog experiments that 
our system is able to maintain reasonable detection capabilities, 
while sampling radar signals that require sampling at a rate of 
about 30MHz at a total rate of lMhz. 

Index Terms — radar, sampling, Compressed Sensing, Sub- 
Nyquist 



I. Introduction 

The classic radar sensing problem treats detection of targets 
moving in space. This is achieved by transmitting RF pulses 
of electro-magnetic energy and sampling the signals caused 
by their reflection. The samples are then processed, in attempt 
to determine the targets' location in space and their velocity. 
Traditional processing methods in both literature and practice 
involve a preliminary stage, referred to as Matched Filtering 
(MF) or Pulse Compression [ 1 ], in which the transmitted pulse 
is correlated with the received signal, in order to increase 
Signal to Noise Ratio (SNR) and improve pulse localization. 
However, such methods are confined to the Nyquist- Shannon 
sampling theorem [2], namely, the received signals must be 
sampled at twice their baseband bandwidth. To achieve high 
spatial resolution and object-distinction capability, modern 
radar systems transmit high instantaneous frequency pulses. 
This implies transmission and detection of wideband signals. 
Consequently, as the demand for resolution increases, classical 
approaches such as MF require sampling at faster rates. This 
implies higher power consumption and vast memory capacity. 

Recently, new approaches |3] to radar processing were 
introduced, which enable sampling and processing at sub- 



Nyquist rates, using analog pre-processing. Related works (4), 
1 5 ] treat ultrasound signal sampling and are readily adapted to 
the radar scenario. These approaches are based on the Finite 
Rate of Innovation (FRI) [6] and Xampling Q frameworks. 
The FRI framework treats sampling and recovery of signals 
characterized by a finite number of degrees of freedom per unit 
time. In many cases such signals can be sampled and recovered 
at a rate proportional to the number of unknowns per time 
interval, which is usually much lower than the Nyquist rate. 
The Xampling philosophy ties together sub-Nyquist sampling 
based on analog pre-processing with techniques of Com- 
pressed Sensing (CS) (8) for recovery. However, as discussed 
in the following sections, these approaches typically require 
sophisticated sampling schemes, which acquire generalized 
measurements of the detected analog signals. Instead, we 
present an analog to digital conversion scheme and a recovery 
algorithm with relaxed constraints, which achieve similar per- 
formance. We further introduce an analog hardware prototype 
which allows sub-Nyquist sampling and reconstruction of 
radar signals, by implementing our proposed scheme. 

Our hardware prototype implements a combintaion of the 
multi-channel topology suggested in and the filtering 
approach presented in [4], additionally taking into account 
the practical challenges they impose. In particular, our board 
consists of 4 channels, each comprising a bandpass crystal 
filter with a random effective carrier frequency. This allows 
to obtain a wide spread of Fourier coefficients of the signal 
in an efficient manner. The proposed recovery algorithm then 
uses these coefficients to recover the delays and amplitudes of 
the radar signal. Using crystal filters, which have extremely 
narrow transition bands, we are able to obtain sufficient 
amount of information from the signal, while substantially 
decreasing the total sampling rate of the channels, as discussed 
in Section [TTT1 Since crystal filters are standard, of-the-shelf 
components, we are confined to adapt our channel design 
to their properties, in order to maximize their efficiency. 
We discuss the challenges this imposes and our method for 
overcoming them in Section Jill 

Simulations as well as real-time experiments of our hard- 
ware prototype prove high detection capabilities for pre- 
integration SNR values over -14dB, with a large reduction 
factor in sampling rate compared to MF. Since practical 
implementation of MF (which makes use of realistic anti- 
aliasing filters) requires oversampling the signal, we have been 
able to achieve a thirty-fold reduction in the sampling rate. 
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That is, our system operates at a total sampling rate of 1MHz 
(20 times less than the signal's Nyquist rate), while an MF 
based system would typically operate at a rate of 30MHz 
(requiring an over-sampling factor of 1.5 in order to use 
practical filters). We are able to achieve this reduction without 
substantially degrading the detection capabilities, provided the 
system operates at high enough SNR. For instance, in our setup 
detection of approximately 90% was achieved at pre-integraion 
SNRs larger than -12 dB. 

To estimate the underlying sparse structure from the Fourier 
coefficients, or the targets in our setting, previous work in 
the FRI framework traditionally employed spectral estimation 
techniques. Here, we use a CS formulation of the recovery 
stage, as suggested in [9]. In CS, the signal is assumed to 
have a sparse representation in a discrete basis. An extensive 
search for the best sparse representation results in combi- 
natorial run-time complexity, which is impractical for real- 
time applications such as radar sensing. Many polynomial-time 
algorithms have been proposed that can be shown to recover 
the true sparse vector under appropriate conditions (8). In our 
simulations, we use the orthogonal matching pursuit (OMP) 
algorithm ITU . 

The remainder of the paper is organized as follows. Sec- 
tion HU establishes the mathematical foundation of our tech- 
niques. We present our radar signal model, link several CS 
concepts to our application and provide justification for our 
chosen parameters. Section [TTT1 compares different analog 
sampling implementations. These systems offer a compromise 
between theoretical requirements and practical hardware con- 
straints. We then elaborate on the design of our hardware. 
Finally, in Section [TV] we present results of MATLAB hard- 
ware simulations, as well as the results of real-time hardware 
experiments. We also describe our realization environment, 
and conclude by discussing the overall system performance. 

II. Xampling of Radar Signals 

A. Methodology 

Radar systems detect targets by transmitting periodic 
streams of pulses and processing their reflections. We model 
the received radar signal as the following stream of pulses: 

M-l L 
n=0 1=1 

where T is the radar's Pulse Repetition Interval (PRI) and M 
is the number of transmitted pulses. This model complies with 
monostatic radars, with non- fluctuating point targets assumed 
to be stationary or moving at very slow velocities. Hence no 
Doppler shift is incorporated into the model. However, future 
work will allow Doppler processing as well. The parameters 
{ai,ti}f =1 correspond to the detected pulses' amplitudes and 
delays respectively, and are proportional to the targets' dis- 
tance from the receiver and their Radar Cross Section (RCS). 
We assume here that the shape of the pulse h(t) and the 
maximal number of pulses L are known, although future 
research may relax this constraint. In particular, knowledge of 
L is used to simplify the stopping criteria of our reconstruction 



algorithm which will be presented in the next part, but is not 
essential. 

Traditional radar systems sample the received signal at the 
Nyquist rate, determined by h (t)'s baseband bandwidth. Our 
goal is to recover x (t) from its samples taken far below this 
rate. It is readily seen that x (t) is completely defined by 
at most 2L unknown parameters, namely a\ and tu within 
any length-T time interval. Hence, in the absence of noise, 
one would expect to be able to accurately recover x (t) from 
only 2L samples per time T 0. If L/T is sufficiently small 
with respect to h(t)'s bandwidth, this implies a significant 
reduction in sampling rate. As shown in (6), 01, H, the 2L 
unknowns defining x(t) may indeed be recovered from only 
2L measurements, corresponding to x(t) 9 s projection onto a 
subset of its Fourier series coefficients. Calculated with respect 
to a single period [0,T), these coefficients are written as 
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where H(uS) is the continuous-time Fourier transform (CTFT) 
of the pulse. Choosing the coefficients such that H (^rfc) is 
nonzero, © can be rewritten as 

X[k] 



Y[k] 
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which is a standard sum-of-exponentials problem. It can be 
shown that this problem has a unique solution given K > 2L 
coefficients 1 6 ] . We discuss the hardware for the acquisition of 
these measurements in Section [Till The rest of this section is 
devoted to reviewing algorithmic approaches for solving ([3]). 

B. Recovery Algorithm 

Many mature techniques for solving ([3]) exist, among which 
are matrix pencil (TU, annihilating filter [6] and many others 
that can be found in [12]. These techniques arise from spectral 
analysis frameworks, and generally require the set of mea- 
surements to form a consecutive subset of the signal's Fourier 
coefficients. 

In it was suggested to use a non-consecutive set of 
Fourier coefficients selected in a distributed manner, as many 
detection systems (namely, ultrasound in (9) and radar in 
our work) benefit from wide frequency aperture. While con- 
secutive coefficients can be obtained using a simple low- 
pass filter, it is shown in [9] that a distributed selection 
results in better recovery and noise robustness. In order to 
remove the constraint of consecutive selection, we choose to 
employ a CS formulation of ([3]) for recovery. We begin by 
quantizing the analog time axis with a resolution step of A, 
thus, approximating ^ as 
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recovery performance is promised, provided that the sens- 
ing matrix satisfies desired properties such as the Restricted 
Isometry Property (RIP). Selecting the frequency samples 
uniformly at random, it is known that if 



Fig. 1: Xampling of radar signals. 

where N = T/A is the number of bins in the PRI. Selecting 
a finite subset of K measurements, namely, k is chosen from 
k = {fci, /u2, . . . , kx} CZ, § may be written as 



y = Ax, 



(5) 



where A is a K x N matrix, formed by taking the set k of 
rows from a N x N DFT matrix, and x is an L-sparse vector 
with non-zero entries at the indices {ni}f =1 . In the context of 
CS, A is known as the sensing matrix. 

Our goal is to find the non-zero entries of x from the 
measurements y. This is a standard CS problem, with A 
being a partial Fourier operator. A solution can be obtained, 
for example, by using the well-known OMP algorithm. The 
algorithm iteratively finds the non-zero entries of x by seeking 
the maximal correlations between y and the columns of A, 
while maintaining an orthogonalization step at the end of 
each iteration in order to improve detectability |[T0lL |8l . 
To improve its performance, we have used a variation of 
the standard OMP algorithm. Specifically, we computed the 
maximal correlation using the pseudo-inverse rather than the 
inner product. Once the nonzero entries are found, the time 
delays are directly calculated and the amplitudes are estimated 
via a standard least-squares technique. A pseudo-code of the 
recovery algorithm is given in Algorithm [T] 

Algorithm 1 OMP Algorithm 

Input: Measurement vector y G Sensing matrix A G 

R XxAr ; Quantization step A; Number of targets L 
Output: Estimated time delays {ti}f =1 \ Corresponding esti- 
mated amplitudes {di}j^ =1 
r <— y; ft <— (j) {Initialization} 
for I — 1 to L do 

p <— A^r {Acquire subspace approximation} 
hi argmax|p n | {Find maximal component} 

n=l,...,N 

ft <— ft U {An J {Augment the estimated delays set} 
A n = 

Pq = 1 — A^A^ {Subspace projection operator} 
r <— P^y {Compute the residual} 
end for 

a = A^A^y {Estimate amplitudes via least- squares} 



[U(fl)'--U(£/)] ]U m (t) = e jT k rnt^ m 



K>CL (log TV) 



(6) 



for some positive constant C, then A obeys the RIP with 
high probability fl3l . In contrast, for consecutive selection 
the RIP is not generally satisfied, unless the cardinality of k is 
significantly increased. However, applying random frequency 
sampling is not practical from a hardware perspective, and 
therefore a rule for selecting a good constellation of frequency 
samples is desired. 

Some practical guidelines for choosing the frequencies are 
suggested in fT4l . The authors consider an equivalent problem 
to ([3]), with the roles of frequency and time interchanged. 
That is, the measurement vector consists of time samples of 
the signal, which obey a spectral line model. Conclusions 
can be easily adapted to our application with minor changes. 
Applying these guidelines, we may formalize the following 
relationship between the support of the recovered signal T, the 
grid resolution A, and the chosen frequency samples {fi}£Li : 
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where c is an empirical constant. In our simulations we set 
c = 20, finding this value to be more robust than that proposed 

in na. 

Condition constraints us to choosing at least two Fourier 
coefficients consecutively. In our application, we choose a 
constellation consisting of four groups of coefficients, where 
the coefficients in each group are consecutive. We refer to this 
choice as multiple band-pass sampling. Note, that © could 
be equally satisfied by simply choosing a single group of K 
consecutive coefficients. As discussed in Section [Till such a 
choice can be implemented by simple hardware, employing 
a single low-pass filtering channel. However, our multiple 
band-pass constellation has the advantage of acquiring the 
measurements over a wider frequency aperture. At the same 
time, it still allows practical hardware implementation (as 
detailed in Section |In|). Referring to ®, by widening the 
frequency aperture we may employ a finer grid resolution 
during the recovery process. Moreover, empirical results show 
that highly distributed frequency samples provide better noise 
robustness |9|. We point out that widening the frequency 
aperture eventually requires increasing the number of samples 
K, otherwise recovery performance may degrade. This trade- 
off, observed in our experiments, is readily seen from ©, 
where the minimal number of samples K gradually increases 
with N, corresponding to the grid resolution. 



C. Frequency Selection 

As mentioned previously, CS based techniques allow flex- 
ibility in choosing the Fourier coefficients. Using OMP, high 



III. Hardware 

We now present practical considerations which guided us 
through the design of our analog board. We begin by examin- 
ing previous proposals for sub-Nyquist schemes and point out 
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the difficulty in their direct implementation. We then present 
our solution. Throughout this section we consider the follow- 
ing system parameters: we are working with a T = 1msec PRI, 
which corresponds to lKHz spacing of Fourier coefficients. 
Our pulse is approximately flat in spectrum, over the extent 
of 10MHz (single-sided band). We assume a maximal number 
of L = 6 targets within the PRI. 

A. Previous Work 

Extracting a consecutive Fourier subset can be done using a 
LPF, followed by sampling at twice its cut-off frequency. The 
DFT of the samples provide the desired Fourier coefficients. 
However, as discussed in Section [HI recovery performance is 
enhanced when using a set of coefficients distributed over a 
larger part of the signal's spectrum. 

For the extraction of arbitrary sets of Fourier coefficients, 
|4] introduces a single-channel, pre-sampling filtering process, 
with the following filter: 

{non-zero, / = ^ , fc G k 
0, J = %,hiK (9) 

arbitrary, otherwise. 

With an arbitrary choice of isolated frequency samples, a 
kernel satisfying © would be characterized by multiple pass- 
bands and extremely high frequency selectivity. In practice, 
these characteristics are difficult to satisfy when designing a 
practical analog filter. 

The selectivity property requires high Q-f actor filters, char- 
acterized by a large attenuation within lKHz, the coefficients' 
spacing. Filtering under this regime requires Q-f actors in 
the order of thousands, which is infeasible even when using 
piezoelectric components. Aside from the problem of achiev- 
ing the high Q-factor, any off-the shelf filter with high Ce- 
faclor will have a long impulse response. This will cause 
instability when processing multiple pulse- streams. Finally, 
there are no off-the-shelf multiple pass-band components. 
These considerations show that implementing this class of 
filters with our application's specification is difficult. 

Another work [ 5 ] suggests the use of multi-channel mixers 
and integrators to directly compute and sample the Fourier 
coefficients at a rate of 1/T in each channel. To obtain 
complex- valued Fourier coefficients, I/Q channels must be 
employed. This doubles the number of channels, and re- 
quires accurate synchronization between the I/Q demodula- 
tors. Furthermore, implementation of multi-channel circuits in 
hardware results in a complicated system, characterized by 
large physical dimensions and high number of components, 
requiring synchronization. 

B. Multi-Channel Crystal Receiver 

Here we present a multi-channel crystal receiver approach 
to obtain Fourier coefficients in a manner that is both practi- 
cal and efficient. This approach makes use of four parallel 
channels which sample distinct bands of the radar signal 
spectrum, as illustrated in Fig. Q] This is achieved through 
filtering the desired band, demodulating it to baseband and 
then sampling it at its Nyquist rate. In this scheme, instead 



of sampling isolated Fourier coefficients, we acquire four 
independent sets of consecutive harmonies. This is a trade- 
off between the theoretical algorithmic requirements, which 
would benefit from a fully distributed selection, and the 
constraints of practical analog filters. In our scenario, the 
proposed approach achieves a 1MHz sampling rate (combining 
all four channels), whereas the Nyquist rate, corresponding to 
h (t)'s bandwidth, is 20MHz. However, the design of systems 
for classical methods such as MF, is constrained to the usage 
of practical filters as well. An anti-aliasing LPF must be used 
prior to sampling, and due to the finite width of its transition 
band, it is not feasible to sample the filter's output at the signal 
Nyquist' s rate without causing aliasing. Modeling a LPF with 
Chebyshev type-I, and allowing a maximal order of 6, we 
were able to achieve a cut-off frequency of 15MHz, which 
requires sampling at 30MHz, a 1.5 oversampling factor with 
respect to the signal's Nyquist rate. Thus we conclude that 
our scheme achieves an even greater reduction compared to 
practical implementations of classical methods. 

In order to maintain a low oversampling factor, we employ 
filters characterized by narrow pass bands. For instance, with 
T = 1msec, a 120KHz pass band corresponds to 120 Fourier 
coefficients. This is ten times the minimal number implied by 
the FRI framework for L = 6 pulses in a noiseless scenario. 
In order to avoid the usage of I/Q channels which, as stated 
before, adds complexity to the system, we demodulate the 
BPF's cut-off frequency in each channel, rather than its central 
frequency, as done in customary design. This requires that the 
filters are characterized by narrow transition bands, in order 
to sufficiently attenuate image frequencies aliasing our chosen 
coefficients. We observed that active filters satisfying our 
narrow pass-band requirement yield unsatisfactory attenuation 
of the image frequencies, degrading the performance of the 
recovery algorithm. We hence chose to use crystal filters, 
whose transition bands are extremely narrow. These filters are 
characterized by a 80KHz pass band. The narrow transition 
bands allow us to achieve low-rate sampling, while extracting a 
sufficient number of Fourier coefficients. Another advantage of 
the crystal filtering implementation is that it results in identical 
receiver channels, which have small phase unbalance and are 
easy to synchronize. 

We performed simulations in order to obtain a good combi- 
nation of the 80KHz pass bands in all four channels: multiple 
constellations of four frequency groups were examined, each 
obtained by randomly choosing four central frequencies and 
then taking narrow sets of consecutive coefficients surrounding 
these frequencies. The following constellation was found to 
yield good performance, based on the evaluation methods 
detailed in Section|TV| 590KHz - 670KHz; 690KHz -770KHz; 
1315KHz - 1395KHz; 1574KHz -1654KHz. 

C Analog Receiver Design 

A full block-diagram of our proposed receiver is depicted 
in Fig. [2l The depicted system is intended to operate in 
baseband frequencies, while the receiving and demodulating 
of RF signals was simulated in a computer environment (see 
Section HVl) . Amplification units were added, in order to 
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Fig. 2: A block diagram of the 4-channel crystal receiver. 



each channel at 250KHz, resulting in a total sampling rate of 
1MHz at all channels combined - a 20-fold reduction relative 
to the pulse's 20MHz Nyquist rate, and a 30-fold reduction 
relative to a practical implementation of MR Once the channel 
output was sampled, its spectrum is easily calculated with the 
FFT algorithm, and the Fourier coefficients are ready as the 
recovery algorithm's input. 

In our design we make a use of 3 filtering stages. The 
first and second stages are located before mixers. When 
using mixers in data paths, special attention must be given 
in order to suppress image frequencies. This is of great 
importance as image frequency coinciding even partially with 
the data Fourier coefficients would result in non-invertible 
distortions. In the first stage modulation to the crystal pass 
band is performed. Since the modulated frequencies and the 
images vary from channel to channel the specification for 
the preliminary LPF must be set to the worst case, that is, 
the left-most image frequency that modulates to our data, 
must be sufficiently suppressed by the LPF. This frequency 
is 56.35MHz. Simulations affirmed that a 40dB attenuation 
is sufficient to obtain good performance, and therefore this is 
the minimal attenuation we require. Table U summarizes the 
specifications for the preliminary filter. 



(a) 



(b) 



(c) 



(d) 



Fig. 3: The dynamics of the signal's spectrum through one 
channel of our 4-channel crystal receiver, (a) Modulation to the 
crystal's frequency, (b) Narrow filtering using a crystal filter. 

(c) The demodulation to baseband and its image frequencies. 

(d) Anti-aliasing filter. 



compensate for the high losses along each channel. Among 
the reasons for these losses, are the four-fold signal power 
reduction at the splitter, the filtering of a 80 KHz slice out of 
10MHz bandwidth signal, and the total four-fold reduction 
caused by the up and down conversions. As crystal filters 
operate only in specific central frequencies, we must perform 
two stages of modulation, in which the desired band is up- 
converted to the pass-band of the filter and then, after filtering, 
is demodulated to baseband. It is then filtered with a 125KHz 
LPF that serves to suppress the modulation image, as well 
as an anti-aliasing filter. The 125KHz cut-off frequency was 
selected in order to acquire over- sampling margin that will 
allow minimal distortion to the Fourier coefficients lying 
in the 25KHz-105KHz band due to image frequencies or 
aliasing. This four-stage process is illustrated in Fig. which 
shows the dynamics of the signal's spectrum through one 
of the receiver's channels. The final LPF allows sampling 



Parameter 


Value 


Pass Frequency 

Maximal Passband Ripple 

Stop Frequency 

Minimal Stopband Attenuation 


2.5MHz 
IdB 

56.35MHz 
40dB 



TABLE I: Preliminary lowpass filter specifications. 

In the next filtering stage, we make use of crystal filters. As 
these are standard, of-the- shelf, devices, the rest of the stages 
must be adapted to their properties, while maximizing the 
channel efficiency (in terms of amount of data acquired and the 
sampling rate at the analog-to-digital converters). Fig. [4| shows 
the magnitude response of a crystal filter, as was measured 
in a network- analyzer, and Table HH details its properties. 
With a data bandwidth of 80KHz at each channel, we are 
able to acquire as many as 320 Fourier coefficients, which 
allows reconstruction of various scenarios including a different 
number of targets, varying distances and a wide RCS range. 
The narrow transition band of the magnitude response, which 
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achieves 60dB attenuation at an offset of 60KHz, allows to 
demodulate the data band to very low frequencies. In practice, 
we demodulate the left-most frequency of the pass-band to 
a frequency of 25KHz, with the right-most demodulated to 
105KHz. Demodulating the data band to the DC frequency 
is undesirable, since in this case image suppression is not 
satisfactory. This is as demodulating the central frequency of 
the crystal filter to DC will require two I/Q channels. While 
this results in a lower sampling rate at each channel separately, 
the total sampling rate will remain approximately the same, 
while the system complexity will increase. 



Parameter 


Value 


Center Frequency 
-3dB Bandwidth 
Maximal Passband Ripple 
Stopband Frequencies 
Minimal Stopband Attenuation 


29MHz 
80KHz 
IdB 

28.94MHz, 29.06MHz 
60dB 



TABLE II: Crystal bandpass filter characteristics. 



Finally, sampling the signal requires that we suppress any 
frequency that might alias to our data. Examining Fig. [3c), 
it is evident that the main contribution to aliasing will occur 
from the part of the crystal's passband energy that was modu- 
lated to very high frequencies, and therefore it is sufficient 
to design an anti-aliasing filter that is efficient enough in 
suppressing that band. However, each channel incorporates 
the usage of amplifiers, which have non-linear regimes. Such 
non-linearities may introduce high-order harmonies of our 
signal, and therefore the anti-aliasing filter must be designed to 
suppress any frequency component which we might not have 
accounted for. Using narrow filters also contributes to noise 
reduction, which improves the performance of the recovery 
algorithm. The desired characteristics of the anti-aliasing LPF 
are detailed in Table Hm With the cut-off frequency being 
125KHz we are able to sample the output signal at a rate of 
250KHz, which yields a total 1MHz rate at all four channels. 



Parameter 


Value 


Pass Frequency 

Maximal Passband Ripple 

Stop Frequency 

Minimal Stopband Attenuation 


105KHz 
IdB 

125KHz 
25dB 



TABLE III: Anti-aliasing lowpass filter specifications. 

Another consideration in the design of the receiver is 
the noise figure (NF) of the system. The NF measures the 
degradation of the SNR from the beginning to the end of the 
channel. In our design, we make use of operational amplifiers 
LMH-6629 as power amplifiers, which have a gain of 14dB 
and an NF of 9dB. The mixers have a conversion loss of 
5dB, while the crystal filters have an insertion loss of 4dB. 
Splitting the power of the signal to four separate channels 
yields attenuation of 6.28dB. The preliminary LPF and the 
anti-aliasing LPF have very small insertion loss, and therefore 
their NF's are negligible. To obtain the total NF of the system, 
we use the Friis equation for the total noise figure of a system 
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Fig. 4: Magnitude response of the Crystal BPF filter as a 
function of frequency. 




Fig. 5: The 4-channel crystal receiver prototype. 



na, 

N F -1 

NF system =NF 1 + Y. rm n r (10) 

where G n and NF n are the gain and noise figure of the 
nth component in the channel, respectively, and N is the 
number of total components. Incorporating the fact that the 
NF of an attenuator is equal to its attenuation level, we are 
able to calculate the total NF of our system, which results 
in NF systern = lldB. This causes a degradation that can be 
compensated in a digital manner, by averaging over several 
periods of the signal, as discussed in Section [iVl 

We present a photo of our 4-channel crystal receiver analog 
board prototype in Fig. In the next section we show real- 
time experiments of our hardware prototype, and compare our 
scheme's performance to the simple LPF channel which filters 
the signal at 500KHz and samples it at an identical rate of 
1MHz, as well as to the traditional MF which operates at 
rates higher than Nyquist. 

IV. Simulations and Experiments 

In this section we present hardware simulations which 
simulate practical filtering, as well as real experiments of our 
hardware prototype that include the transmission of analog 
signals. 

A. Simulations 

The first evaluation step we introduce is simulating our 
hardware and recovery method via MATLAB. We examined 
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our technique's success in recovering multiple realizations of 
the signal defined in (Q]), after its corruption by noise. Each 
realization comprised L = 6 pulses. The pulse and its spec- 
trum are depicted in Fig. [6] Time delays and amplitudes were 
drawn uniformly, at random, within the intervals [0,1msec) and 
[0.5,1.5] respectively. The signals were generated digitally, at 
a rate of 0.6GHz, which is much higher than the Nyquist rate 
corresponding to h(t) 9 s 10MHz bandwidth. The signals were 
corrupted by zero-mean white Gaussian noise, with variance 
a 2 determined such that the SNR, defined with respect to the 
weakest target as 



SNR 



1 • , ,2 

— mm \ai\ 



(11) 



maintains a predefined value. 

To measure the system's recovery performance and to 
compare it to other solutions, we define the Hit-Rate and the 
RMSE metrics in the following manner: 



Hit-Rate = - \{U\ \i t - t t \ < e th ,l = 1,...,L}\ 



RMSE 



1 



V 



L • Hit-Rate 



E $-*Y 



(12) 
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where {ii}f =1 are the estimated time delays and e t h is the 
tolerance factor, determined by application (we chose e t h = 
150ns, which is three times the signal's Nyquist period). The 
RMSE is calculated only with respect to the estimations found 
within the tolerance interval. In order to obtain statistically 
stable results, each experiment was repeated 500 times. 

We evaluated the performance of four different estimation 
methods. The first one was our 4-channel crystal receiver, 
where OMP was used for recovering the unknown signal 
parameters from its subset of Fourier series coefficients. The 
next two methods were based on a single LPF channel which 
operates in a total sampling rate of 1MHz, and obtains a group 
of consecutive Fourier harmonics. Although the sampling rate 
at this single channel is 1MHz, we are only able to extract 
400 coefficients due to non-ideal transition bands. In one 
of the methods, a spectral analysis method, known as the 
matrix pencil (MP) algorithm ifTTl was used, while in the 
second we used again the OMP algorithm. Note that all three 
of the schemes operated under the same overall sampling 
rate. Finally, we compared the former approaches with the 
traditional MF method, which requires sampling at a rate 1.5 
times higher the Nyquist rate (30MHz). 

All the filters used in our simulations were modeled via 
MATLAB's filter design and analysis tool, and were based 
on the Chebyshev type-I filter. The electronic filters in the 
different schemes (such as the anti-aliasing LPF's) were de- 
signed so that their order would not increase above 6. In 
contrast, the crystal filter was modeled with a high order 
BPF, in order to mimmick its narrow transition bands. As 
part of the digital procedure, we integrated over reflections 
of 500 consecutive pulses. This was done in order to reduce 
noise power and improve detectability. While this procedure 
constrains us to slow moving targets, our technique remains 



(a) 



(b) 



Fig. 6: The temporal shape of the pulse (a) and its Fourier 
spectrum (b). 



applicable to various scenarios, such as naval target detection, 
vehicle and human tracking, etc. 

Fig. |7] depicts the performance of the four aforementioned 
methods, as a function of the pre-integrated SNR, tested for 
multiple realizations of our radar signal. Note that two of the 
methods make use of the same single LPF channel. Examining 
the results we infer that our 4-channel crystal receiver yields 
much better performance than LPF based methods at lower 
SNR values, an outcome that makes sense as in noisy realiza- 
tions the aperture is critical. As the SNR increases though, the 
aperture plays a less significant role, and the total number of 
samples determines the reconstruction performance. It is not 
surprising that the less noisy the samples are the better the 
performance is obtained by MP. The latter does not quantize 
the time axis and is proved to reach numerical precision for 
high enough SNR values IfTTl . For complete comparison, the 
MF curve is plotted. It is quite clear that it provides better 
performance than our approach, in particular for lower SNR 
values, but as the SNR improves the performance gap drops. 
We also infer that our 4-channel crystal receiver provides lower 
RMSE than the LPF channel based methods, in the entire 
SNR range, and it gets closer to the MF's performance as the 
SNR increases. Therefore, we can conclude that our system 
offers the best recovery performance apart from MF (which 
requires fast sampling rate by Nyquist) in the SNR range that 
is effective for realistic radar scenarios. 

B. Experimental Environment 

To evaluate the board we make use of NI PXI equipment 
for both system synchronization and signal sources. The entire 
component ensemble, wrapped in the NI chassis, is depicted in 
Fig. [8] The components we make use of are listed in Table HVl 

In order to achieve system stability, and to obtain fine 
recovery performance, synchronizing the system is essential. 
Jitters and drifts between clocks, might skew the time instances 
in a manner that distorts the calculated Fourier coefficients. 
We therefore must ensure that our entire component ensemble, 
including the AWG, LO signals and ADCs, are triggered by 
the same clock. This promises that phase unbalance between 
devices is brought to a minimum. 

The advantage of the NI PXI chassis is that it can synchro- 
nize several devices to one base clock and distribute a trigger 
signal with skew less than Ins. The ability to synchronize 
all devices to one clock keeps the jitter between signals low 
so that we can assume it does not cause a change in phase 
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Fig. 7: A comparison of four detection methods, (a) Hit Rate 
(b) RMSE. 



during the operation. Therefore, the whole system is stable 
and consistent. 



NI Part No. 


Description 


NI PXIe-1075 
NI PXIe-8133 
NI PXIe-5451 
NI PXI-5690 
NI PXI-4130 
NI PXI-6123 

NI PXIe-6672 

NI 5781 

NI PXI-7965R 


18-Slot 3U PXI Express Chassis 
Core i7-820QM 1.73GHz Controller 
400MSamp/s Arbitrary Waveform Generator (AWG) 
2 Channel RF Preamplifier from 500KHz to 3.0GHz 
Power SMU 

16-Bit, 500KSamp/s/ch, Simultaneous Sampling Multi- 
function ADC 

Timing and Multichassis Synchronization Unit for PXI 
Express 

Baseband Transceiver for NI FlexRIO 
NI FlexRIO FPGA Module for PXI 



TABLE IV: NI PXI devices. 



C. Experiments 

The experimental process consists of the following steps. 
We begin by using the AWR software, which gives us the 
ability to examine a large variety of scenarios, comprized 
of different targets, distances and RCS values. It is able 



to simulate the complete radar scenario, including the pulse 
transmission and accurate power loss due to wave propagation 
in a realistic medium. It also takes into account the reflec- 
tions from the detected objects, which are proportional to 
the targets' RCS. The AWR software also contains a model 
of a realistic RF receiver, which performs signal processing 
in high frequencies. The demodulation of the signal to IF 
frequencies is simulated, and the output is saved to a file. 
Next, the simulation result is loaded to the AWG module, 
which produces an analog signal. This signal is amplified using 
the NI 5690 LNA and then routed to our 4-channel crystal 
receiver. The receiver is fed by 5 LO, of which 4 modulate 
the desired frequency band of each channel individually to 
the crystal pass band, and a global one that modulates the 
latter to a low frequncy band, before sampling. The LOs 
are created using 3 NI 5781 baseband transceivers, acting as 
trigger based signal generators with a constant and known 
phase, controlled by NI Flex Rio FPGAs. The AWG also 
triggers the ADC to sample 250 samples in each sampling 
cycle, per channel. These samples are fed into the chassis' 
controller and a MATLAB code is launched that runs the 
reconstruction algorithm. Our system contains a fully detailed 
interface implemented in the Lab View environment, which 
allows simple activation of the process. A screenshot of the 
interface is depicted in Fig. 

Performing the previously described processes, we tested 
various scenarios on the board - that is, a variety of targets, 
distances and RCS - and examined the reconstruction quality. 
We tested our board for its limitations and have obtained the 
following results. The minimal voltage for a reconstructible 
pulse, without PRI integration, is approximately 4mV, which 
corresponds to a power of P m in = 320nW = — 34.94dBm. As 
the board contains amplification units, we must also consider 
the maximal allowable input power. Above this power, the 
amplifiers reach their non-linear regimes and introduce non- 
invertible distortions to the input signal. It was measured that 
above a voltage of 92mV, which corresponds to a power of 
Pmax = 169.28/iW = -7.7139dBm, such effects begin to 
take place. We thus deduce that our system has a dynamic 
range (DR) of 



DR 



system 



27.23dB. 



(13) 



The DR of the analog system can be digitally enhanced using 
pulse integration. For example, integrating over 500 pulses, 
increases the DR by a factor of 101og 10 (500) = 27dB. 

Moreover, we have tested the limitations for multiple- 
pulse scenario. In cases where the OMP algorithm fails to 
obtain time delays with a small enough error, calculating 
the residual might yield unsatisfactory attenuation of the part 
of the measurements corresponding to the located pulse. For 
pulses energetic enough this means that their residual might 
still be stronger than the remaining pulses. We have measured 
that a scenario of two pulses of which voltage ratio is above 
12, causes the algorithm to select the stronger pulse twice in 
consecutive iterations, thus failing to identify the weaker pulse. 

In Fig. [TOl we show a realization of a radar signal with 
L = 4 targets and its reconstruction using our 4-channel 
crystal receiver. In order to achieve a low SNR scenario, we 
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Fig. 8: The NI Chassis. 




Fig. 9: The Lab View experimental interface. From left to 
right: the signal's spectrum, the frequency response of each 
channel, the 4 signals detectd in each channel, at the top- the 
reconstructed signal, at the bottom- the transmitted signal. 



have reduced the pulses' power to approximately 5dB above 
the system's noise level. The input noise comes mainly from 
the AWG's thermal noise, and has a power of approximately 
— 75dBm. In order to enhance the SNR, an averaging over 100 
periods of the sampled signal was performed prior to the FFT 
calculation. It can be seen that the time delays were recovered 
with high precision. However, there is a slight error in the 
estimation of the amplitudes. 

V. Conclusions 

We presented the first sub-Nyquist radar prototype based 
on the Xampling methodology. We have been able to reduce 
the total sampling rate by a factor of approximately 30 while 
maintaining reasonable detection ability. 

Our experiments with the hardware have further proven 
that sub-Nyquist sampling of radar signals is possible without 
losing much of the recovery performance. 

In the future, theoretical research might enrich the mathe- 
matical model, with the addition of Doppler effect and relaxing 
the demand that the shape of the pulses is known. Different 
algorithms and methods might be considered instead of OMP, 
so that noise robustness and spatial resolution might increase. 
However, we believe that the actual implementation of the 
sub-Nyquist radar prototype, is an important and crucial step 
towards the incorporation of Xampling and FRI frameworks 
into real communications and signal-processing systems. 
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Fig. 10: (a) A typical radar signal corrupted by AWGN and 
(b) the noise free signal and its recovery. Pre-integration SNR 
is -lOdB. 
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